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Using D* + mesons exclusively reconstructed in the DELPHI detector at LEP, orbital and radial excitations of non strange charmed 
mesons are studied. The multiplicities of the two narrow Z)J and D^ orbital excitations are measured in Z° — » cc and Z° —> bb 
decays. Preliminary results are obtained on the production of broad D** states, using B meson semi-leptonic decays. A narrow 
signal of 66 ± 14 events is observed in the (D* + 7r + 7r~) final state, interpreted as the first evidence of the predicted D* radial 
excitation. 



1 Introduction 

For mesons containing heavy and light quarks (Qq), and 
in the limit where the heavy quark mass is much larger 
than the typical QCD scale, the spin sq of the heavy 
quark decouples from other degrees of freedom. Thus, 
for strong decays, the total (spin+orbital) angular mo- 
mentum j q = s~^+ L of the light component is conserved. 
This heavy quark symmetry, together with quark poten- 
tial models used for lower mass mesons, allows the masses 
and decay widths of heavy mesons to be predicted^. 
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Figure 1: Spectroscopy of non-strange D mesons. The shaded areas 
show predicted widths for these states. For clarity the expected D± 
and decays involving a p meson or tttt pairs are not shown. 



The present knowledge of charmed meson spec- 
troscopy is summarized in Figure 0. The well established 
D and D* mesonsB correspond to the two degenerate lev- 



els of the (L=0, j q = 1/2) state. The two (L=l, j q = 3/2) 
states have been clearly observed 13, because they have 
narrow decay widths of about 20 MeV/c 2 . The mea- 
sured masses of the D?(2420) and D2°(2460) agree within 
20 MeV/c 2 with the prediction of the models. Section 3 
presents a measurement of their production rate in cc 
and bb jets. 

The (L=l, j q = 1/2) states decay through a S wave 
and are expected to have large decay widths. Up to now, 
they have not been observed directly, but their total pro- 
duction rate is measured using B meson semi-leptonic 
decays (section 4). 

In addition to these orbital excitations, ra- 
dial excitations of heavy mesons are foreseen. 
The D and D* are expected to have masses of 
2.58 GeV/c 2 and 2.64 GeV/c 2 respectively, with j& 
10-25 MeV/c 2 uncertainty on the mass predictions 13. 
They are expected to decay, in S wave, into D^tt-k. 
Section 5 presents the first evidence for the D* meson, 
observed in the decay mode (D*mr). 

2 D** and D* reconstruction 

DELPHI @ is a multipurpose LEP detector, with spe- 
cial emphasis on precise vertex and charged tracks mo- 
mentum reconstruction, and particle identification. The 
micro-vertex detector provides 3 and 2 Z hits per 
track, with intrinsic resolutions of 7.6 and 9 /im. For 
muons of 45 GeV/c momentum, a resolution of a(p)/p of 
± 3% is obtained, and the precision of the track extrap- 
olation to the beam collision point is 26 ± 2 /an. Kaon 
and pion identification is performed using a Ring Imag- 
ing CHerenkov detector, and the ionisation loss in the 
TPC, which is the main tracking device. A total of 3.4 
million hadronic events is obtained from the 1992-1995 
data, at center-of-mass energies close to the Z° mass. 

2.1 D* reconstruction 

All the decay channels considered here involve the 
D*+ -> D°ir+ decay, followed by D° -> {R--K+) or 
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Figure 2: A"7r and i\37r invariant mass distributions in (a) the 
D*+ -* (ftT-7r+)7r+ and (b) the D*+ (ft:-7r+7r-7r+)7r+ decay 
channels. 



To reconstruct the D° decay final state, all (K~ir + ) 
and (K~ir + ir~ir + ) combinations are tried to fit a sec- 
ondary vertex in space. Kinematieal and track selec- 
tion cuts are described in detail in la. Kaon candidates 
are considered if they have a momentum larger than 
1 GeV/c and, in the K3ir channel, a loose kaon iden- 
tification is required. The D° momentum and invari- 
ant mass are computed from the momenta of the decay 
products. Then, all charged particles with momentum 
between 0.4 GeV/c and 4.5 GeV/c and charge opposite 
to that of the kaon candidate are used as pion candidates 
for the D*+ -> D°n+ decay. In the Kn (K3n) channel, 
events are selected if the mass difference (Mktt-k, —Mk-k) 
(resp. (Mic 37r7r , -M K3 „)) is within ± 2 (± 1 ) MeV/c 2 of 
the nominal value (Mr>* — M D o). The D* candidates 
must have an energy fraction Xe(D*) — E(D*)/Eb eam 
greater than 0.25. Figure ^| shows the distribution of the 
M(Ktt) and M(K3ir) invariant masses for the selected 
events. The fitted D° masses and widths are 1868 ± 1 
(1869 ± 1) and 19 ± 1 (12 ± 2) MeV/c 2 . The recon- 
structed D° mass is required to lie within ± 40 (± 30) 
MeV/c 2 of the nominal D° mass: 4661 ± 88 (2164 ± 65) 
D* candidates arc selected in the Ktt ( K3ir) channels. 
The selection efficiency is estimated, using the simula- 
tion, to be 21% (8%). 

2.2 D\, D^ and D* reconstuction 

Similar selection criteria and vertex reconstruction are 
used to reconstruct narrow orbitally and radially excited 
states. 



a Throughout this paper, charge-conjugate states are always im- 
plied, and the it from D* decay is denoted as 7r*. 



In the case of and DZr decaying into D* + ir~ , a 
pion with a charge opposite that of the D* + is added, and 
the D°7r+7r~ vertex is fitted. All combinations are tried, 
provided the pion candidate has a momentum larger than 
1.0 (1.5) GeV/c in the Ktt (K3w) channel. The recon- 
struction efficiency is 14% (6%) in the Ktt ( KZtt) chan- 
nels. 

In the case of D* decaying into D* + 7r + 7r~, all pairs of 
oppositely charged pions are used to fit a D°tt + tt~ vertex. 
The pion candidates are required to have a momentum 
larger than 0.6(1.0) GeV/c, and those compatible with 
a kaon according to particle identification are rejected. 
For a signal of mass 2640 MeV/c 2 , the reconstruction 
efficiency is 4% (2%) in the Ktt ( K3ir) channels. 

In both cases, the precision on the invariant 
mass reconstruction is improved by correcting for a 
4 MeV/c 2 shift observed in the D° mass, by using: 



M{D*tt) 
M{D*tttt) : 



M, 



(-D°7T*7r7r) 



M, 



(£>°7T») 

(■D°7r„) 



m D * 
mo* 



(1) 



where m/j. is the nominal D* + mass. The simulation 
predicts a resolution of about 6 MeV/c 2 on the mass re- 
construction, for both radial and orbital excitations. 

2.3 Selection of bb and cc samples 

Due to the relatively long lifetimes of charmed and bot- 
tomed particles, heavy flavour events are characterized 
by the presence of secondary vertices. The probabil- 
ity V that all tracks detected in the event come from 
the primary vertex is small: for bb events, a purity of 
90% is archieved, with an efficiency of 60%, by requiring 
V < 10~ 2 . 

Charmed mesons from Z° —* bb events are distin- 
guished from those in cc events by considering both their 
energy and lifetime informations. Bottom quarks frag- 
ment into a B hadron, which subsequently decays into 
a D* + meson, whereas in cc events charmed mesons are 
directly produced in the fragmentation process. This dif- 
ference in the hadronization leads to a smaller energy 
fraction of Xe(D*) for bb events. Also, due to the b 
quark lifetime, the apparent flight of the D° meson is 
greater than the true decay length. Its measured proper 
time distribution is larger than the mean B meson life- 
time, 1.6 ps, compared to a true D° lifetime of 0.4 ps. 

By combining these variables, bb and cc samples are 
selected, with high purities: 92 % for bb , 89 % for cc. 

In the bb sample, the combinatorial background is 
higher, but is reduced by 50% using the kaon identifica- 
tion, and also by asking that the impact parameter of 
the additional pion is positive, i.e. that the intersection 
of the pion and D* directions is on the same side of the 
primary vertex as the D* vertex. As a consequence, the 
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ratio of efficiencies e{D* ir) / e{D*) is 52% for bb , com- 
pared to 62% for cc . 



3 Study of narrow orbital excitations 

Figure [| shows the M(D*tt) invariant mass distribu- 
tion obtained for the sum of the bb and cc samples B. 
A clear excess of (D* + tt~) pairs is observed between 
2.4 and 2.5 GeV/c 2 , corresponding to the two overlap- 
ping contributions of the D\ and Df. They are fitted 
by two Breit-Wigner functions, .whose widths are fixed 
to the measured world average convoluted with the 
experimental resolution. A total signal of (361 ± 58) 
D° + Df events is fitted, out of which (65 ± 10) % 
is assigned to D®. The masses are left free in the fit, 
and the result is M D o = 2425 ± 3(stat) MeV/c 2 and 
M D ,o = 2461 ± 6{stat) MeV/c 2 , i.e. consistent with the 

world averagesi. The helicity distributions are consistent 
with the production of a J p = 1 + and J p = 2 + states. 

Figure |] shows the same mass distribution, but for 
the bb and cc samples separately. The same fit is per- 
formed, but both D® and Df masses and widths are 
fixed to the world average. The result of the fit is 
(97 ± 26) Dl and (69 ± 27) Df in the bb sample, 
(141 ± 26) Dl and (104 ± 26) Df in the cc sample. 
In order to measure the D® and Df production rates, 
these results are unfolded from the reconstruction effi- 
ciencies, signal purities, D® and Df decay widths into 
D* + ir~. The errors quoted below are statistical only. 
Systematic errors are still under study, but smaller than 
the statistical errors. 

The cc sample provides direct information on the 
charm fragmentation. Results are: 



/(c -> D\) = 1.9 ± 0.4 {stat) % 
f(c Df ) = 4.3 ± 1.3 {stat) % 



(2) 



Both results are in agreement with previous LEP and 
CLEO measurements. For the the result is also 
in agreement with theoretical calculations 13, which pre- 
dicted 1.7%. For the Df, the result is high compared to 
the expectation (2.4 %), but has large errors. A more pre- 
cise measurement would need to use the Df — > D + ir~ 
channel, which is forbidden for the D®. 
For the bb sample, results are: 



Df 



2.2 ± 0.6 {stat) % 
) = 4.8 ± 2.0 {stat) % 



(3) 
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Figure 3: Invariant mass distributions (D*+7r — ) (dots) and 
(D* + 7r + ) (hatched histogram). The mass computation and the 
fit are explained in the text. 
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This shows that the charm fragmentation properties are 
similar, although in a different environment. 



Figure 4: Invariant mass distributions (D*+7r — ) (dots) and 
(D* + 7r + ) (hatched histogram) for the separate bb and cc samples. 
The mass computation and the fit are explained in the text. 
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4 Study of broad orbital excitations in B meson 
semileptonic decays 

In B meson semileptonic decays, only 60% to 70% of 
the final states are described by DiDi and D*£vg. The 
remaining contribution is attributed to D**. The total 
production rate, including broad orbital excitations, can 
be measured using the impact parameter of the pion, 
denoted 7r*», emitted in the decay chain B — > D**X — > 

(D*7T**)X'. 

Events are selected if a lepton with momentum larger 
than 3 GeV/c 2 is identified, and if its transverse momen- 
tum relative to the D* + is larger than 0.5 GeV/c 2 . The 
kaon candidates in the D° decay must have the same 
charge as the lepton. 459 ±25 (288 ± 19) events are se- 
lected in the Ktt(K3tt) channel. All remaining tracks, of 
charge opposite to that of the D* + , are 7r*» candidates. 
Kinematical and selection cuts are described in B. The 
background due to fake D* + associated to a true lep- 
ton £~ is subtracted by using events in the tail of the 
D* invariant mass distribution. The contribution of true 
D* + associated to a fake lepton is subtracted using D*£ 
pairs with the wrong sign combination. The remaining 
111 ± 16 events are due to true b semileptonic decays 
into D* + £~X final state, associated with a 7r»* candi- 
date either from D** decay, or from jet fragmentation. 
The shape of the two contributions are shown in figure ||. 
They are used to fit the 7r*„ impact parameter distribu- 
tion shown in figure ^. From the result of the fit, the 
following branching ratio is obtained: 



BR(B- -> {D*+tt-)£-9X) 
1.15 ± 0.17(stat) ± O.U(syst) % 



(4) 
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Figure 5: Impact parameter relative to the primary interaction 
vertex in simulated B semileptonic decays for a) 7r** from D** 
decay and b) charged particles from jet fragmentation. 
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Figure 6: Impact parameter relative to the primary interaction 
vertex in real data for a) right charge 7rJ"„ and b) wrong charge 
7T«* candidates. The hatched area is the contribution from jet 
fragmentation. 



This result significantly improves a previous 
DELPHI measurement, and is in agreement with other 
LEP measurements. 

The (D* + 7r»*) invariant mass is also reconstructed 
and used to fit, in the way described in the previous 
section, the and D^ narrow resonances. A signal of 
26.7 ± 8.2 D\ is fitted, and the corresponding production 
rate is: 



BR(B- -> D\£~vX) 
0.72 ± 0.22(stat) ± Q.13(syst) % 



(5) 



For the state, 14.8 ± 7.7 events are fitted, i.e. a 
signal significance smaller than 2 a. More data would be 
necessary to estimate the corresponding branching ratio. 



5 Evidence for a narrow radial excitation 

Figure shows the invariant mass distribution obtained 
when two pions of opposite charges are added to the D* 
candidate, and using the sum of the two bb and cc sam- 
ples. An excess of 66 ± 14 (stat) events is observed in the 
(D* + ir + ir~) combination. The signal is fitted by a Gaus- 
sian distribution of free parameters: the x 2 P er degree of 
freedom is 60/59, and would be 91/62 if the Gaussian 
was removed. About (57 ± 10)% of the signal is selected 
in the cc sample. 

The fitted mass is 2637 ± 2 (stat) ± 6 (syst) MeV/c 2 . 
It is thus consistent with the predictions for the D* ra- 
dial excitation I: 2640 MeV/c 2 . Other L=2 states are 
predicted, but with masses higher by at least 50 MeV/c 2 . 
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The width of the fitted Gaussian is 7 ± 2 MeV/c 2 , i.e. 
compatible with the detector resolution. Therefore, only 
an upper limit is derived: the full width of the signal is 
smaller than 15 MeV/c 2 at 95 % C.L. There is no natural 
explanation of such a small value, neither for the D* nor 
for higher orbital excitationsla. 

Various checks were performed. Varying the back- 
ground shape and the kinematical cuts has no effect 
within statistics. No peculiar double counting was no- 
ticed, and the signal is stable when the tt* is added to 
the D°tt + tt~ tracks in the vertex fit. As explained above, 
mass shifts are studied using D\ and D^ narrow states, 
and a conservative systematic error of 6 MeV/c 2 is at- 
tached to the mass measurement. 

The production rate of this signal can be compared 
with that of the D\ and D^ narrow states: 



The measured multiplicities in bb events are consistent 
with the ones in cc events, both for the D\ and D^ ■ 

The total D** production rate, involving a D* + in the 
final state, is measured in B meson semileptonic decays. 

A narrow signal is observed in the (D* + ir + ir~) fi- 
nal state, at the mass M = 2637 ± 2 (stat) ± 6 
(syst) MeV/c 2 , interpreted as the first evidence of the 
predicted D* meson. 
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<N D ,,>xBr(D* ->D*7r + iO 

E J=1 ,2 <Ar D (;) >xSr(D( /'^ D * 7r) 
0.49 ± 0.18(stat) ± O.lO(syst) 



(6) 



Most of the systematic uncertainties cancel in this 
ratio. The quoted systematics is due to the Monte-Carlo 
statistics, and to the uncertainties on widths and on the 
kaon rejection. This result is compatible, within its large 
errors, with the value obtained using the thermodynam- 
ical models already mentioned for orbital states Bfl. 
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Figure 7: Invariant mass distributions (D*+7r+7r _ ) (dots) and 
(D* + ir~ tt~ ) (hatched histogram). The mass computation and the 
fit are explained in the text. 



6 Conclusion 

Using about 7000 exclusively reconstructed D* mesons, 
the D® and D^ multiplicities are measured in cc events, 
and found to be consistent with theoretical calculations. 
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